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We investigate the effect of hole doping on the strong-coupling Hubbard model at half-filling 
in spatial dimensions D > 1. We start with an antiferromagnetic mean-field description of the 
, insulating state, and show that doping creates solitons in the antiferromagnetic background. In one 

■ dimension, the soliton is topological, spinless, and decoupled from the background antiferromagnetic 

fluctuations at low energies. In two dimensions and above, the soliton is non-topological, has spin 
quantum number 1/2, and is strongly coupled to the antiferromagnetic fluctuations. We derive 
' the effective action governing the quasiparticle motion, study the properties of a single carrier, and 

I comment on a possible description at finite concentration. 

CSl ■ PACS Numbers: 71.10.Fd, 74.20.Mn, 71.27.+a 



^ • I. INTRODUCTION 

O 

o\ . 

<^ , Nearly sixty years ago, Peierls first suggested that correlation effects are the cause of the insulating behavior in 
C\| ■ nickel oxide.B Today, we know that NiO is an example of a "Mott insulator" - an insulator that would be metallic by 
valence counting. The parent compounds of high temperature superconductors, e.g., La2Cu04, are well described as 
^n , Mott insulators, and have stimulated a renewed interest in the subject. 

Mott insulators are generally antiferromagnetic, a feature which makes a mean-field spin-density-wave (SDW) 
description possibleH. Although an antiferromagnetic insulator within the SDW description resembles a band insulator, 
the difference between the two types of insulators is most pronounced when charge carriers are added. Doping a band 
5^ insulator (say with p-type dopants) moves the Fermi energy from midgap into the valence band. The hole "pocket" in 
the valence band is responsible for the conductivity. In a Mott insulator, however, doping creates solitons: localized 
^ regions of suppressed antiferromagnetic order in which the hole is confined. The energies of the soliton states lie inside 
O ' the SDW gap and pin the chemical potential. The soliton propagation is responsible for the conductivity. 
J-^ ' The nature of the solitons in Mott insulators depends on therSaatial dimensionality. In one dimension (ID), the 
^ , solitons are topological - they are the anti-phase domain wallE'u They have spin quantum number zero and are 
' decoupled from the spin fkictuations at low energies. As a result of the decoupling, the quasiparticle bandwidth is 
rS only weakly renormalized.Q In D > 2, the solitons are non-topological ~ they are "bag"-li^je|-|Objects in which the 
' antiferromagnetic order is suppressed consistent with the symmetry of the underlying lattice.&Ll They carry spin and 
interact strongly with the low-energy spin fluctuations, significantly reducing the quasiparticle bandwidth. 

Both numerical evidence and the behavior of the optical conductivity in the high temperature superconductors are 
suggestive of the existence of such quasiparticles in two dimensions- Numerical solutions for the single particle Greenia 
function indicate a single pole and an incoherent backgroundH" Il3 In the high Tc oxides, the optical conductivitytil 
contains two components, a Drude peak and a broad mid-infrared absorption, which are related to the added charge 
carriers (their spectral weight satisfies the sum rule) and can not be well fit over the entire region assuming conventional 
quasiparticles. Furthermore, the two components have very different temperature dependence. While the Drude peak 
is sensitive to temperature (e.g. the width of the Drude peak is « 2A:bT), the mid-infrared part is not. In fact, 
the latter survives even in the superconducting state. This type of optical conductivity is, however, what would be 
expected from solitonic quasiparticles whose conduction is responsible for the Drude peak and whose internal structure 
generates the mid-infcired absorption. Such behavior is seen in other doped Mott insulators,ll3 and systems such as 
doped PolyacetylenctJ in which solitons are known to exist. 

In this paper, we solve the mean-field equations self-consistently to determine the internal structure of the soliton 
and then examine the motion of the soliton in the antiferromagnetic background. In the mean-field treatment described 
below, the directional fluctuations of the order parameter are turned off. Consequently, the mean-fleld solution exhibits 
long-range order which breaks the lattice translational symmetry. At low energy and long wavelength the important 
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dynamical degrees of freedom are a) the motion of the charged sohtons, and b) the directional fluctuation of the 
antiferromagnetic order parameter. These quantum fluctuations restore the symmetry broken at the mean-field level. 
To obtain an effective action governing these important dynamical processes, we force the antiferromagnetic order 
parameter to have both solitons (with time-dependent locations) and space-time-dependent directional twists. We 
then eliminate the electron degrees of freedom corresponding to particle-hole pair excitations. 

In one dimeasion, the soliton in the Hubbard model is intimately related to the structure soliton in doped 
Polyacetylencim In Polyacetylenc, the soliton arises because the ground state of the undopcd material breaks a 
two-fold symmetry of the total electron-phonon Hamiltonian. For the undoped (i.e. half-filled) Hubbard model in ID, 
there is no known broken symmetry. The spin-spin correlation decays algebraically.li3 From this point of view it is 
surprising that there exist finite-size solitons upon doping. The key to the solution of the puzzle is the realization of 
the very different roles played by the longitudinal (magnitude) and transverse (direction) components of the antiferro- 
magnetic order parameter. A useful picture for the algebraic order in the ID Hubbard model is the following: Imagine 
locally (in the sense of both space and time) the order parameter has a well-developed magnitude. However, on a 
much larger space-time scale, the direction of these local order parameters fluctuates, which destroys the long-range 
order. In other words, the algebraic antiferromagnetic correlation is due to the power-law decay of the directional 
correlation between local order parameters, not due to-the absence of local moments. In the case of high Tc oxides 
(2D system), this picture is supported by experiments£3 Thus the mean- field longitudinal spin order plays precisely 
the same role as the Peierls order does in Polyacetylenc. 

The soliton in ID is a 7r-phase-shift domain wall between the two different mean-field ground states. In the absence 
of doping the soliton is unstable. In fact once we allow the spin direction to fluctuate, the soliton decays into an 
infinite-width spin twist. Finite-size solitons are stabilized by doping. As will be shown later, after accommodating a 
hole, a charge soliton remains sharp even when the directional fluctuations are allowed. The resulting action governing 
the charge and spin dynamics are given by Eqs. ( p^ ) and (|3^). Here we note that due to the fact that the charged 
solitons are spinless, the two degrees of freedom are decoupled. 

In 2D a domain wall is an extended object. Thus one may expect that it is more difficult for doping to stabilize 
it. For one hole in a half- filled Hubbard model with sufficiently big U/t [t is the hopping amplitude and U is the 
onsite repulsion strength), we find that this expectation is indeed true. In this case we find that the doped charge is 
accommodated by a localized non-topological soliton (or a spin bag). This soliton is the finite U generalization of a 
hole in the t-J model: It has spin quantum number 1/2 and its spin direction is constrained to be opposite to the 
magnetic moment in its absence. Here due to the finiteness of U , a small amount of double occupation exists, and the 
size of the soliton depends on U/t. The effective action describing soliton motion and antiferromagnetic fluctuations 
is the t-J model. Due to the coupling to the antiferromagnetic fluctuations, the soliton will be further dressed to 
form a quasiparticle. 

The quasiparticle formation can be-siHip|ly understood in the string picture by considering the soliton propagation 
in the antiferromagnetic background.LHlj'E2l As the soliton moves, it leaves behind a string of overturned spins which 
cost antiferromagnetic energy. This energy cost inhibits the motion of the soliton and mimics the effects of a linearly 
confining potential. The bound state of the soliton in this potential is the quasiparticle E3 It is important to point out 
that spin ffipping processes in principle can erase the string of frustrated bonds behind the soliton and enable it to 
escape the potential. If that happens the quasiparticle is no longer well defined. Exact numerical studies of one hole in 
an antiferromagnet, however, suggest that this does not happen and the quasiparticle picture remains valid. For the 
2D Hubbard model with t ancLJi-chipseji in the range relevant to the high- Tc superconductors, the quasiparticle size 
is about 3 ~ 5 lattice spacings.ljii30~Ea We therefore expect that the qualitative properties of the quasiparticle, such 
as its quantum numbers, are not changed by the omitted spin flipping processes. However, all quantitative aspects 
are subject to renormalization. In addition to their effect on the internal structure of the quasiparticle, quantum spin 
fiuctuations have another important effect - they generate coherent propagation of the quasiparticle as a whole. The 
result is that the quasiparticle primarily hops between neighboring sites within the same sublattice. In the presence 
of slow time-dependent twisting of the antiferromagnetic order, the quasiparticle spin must follow the direction of 
the instantaneous local moment. Thus we arrive at the following "surfing model": Imagine a slowly fiuctuating 
membrane whose normal directions define the local directions of the magnetic moments, then the quasiparticles surf 
on the membrane but have to keep their spins locally perpendicular to the membrane. This is illustrated in Fig. ^ It 
is knownE3n23 that when a quantum mechanical particle moves on a fluctuating curved surface under the constraint 
that its spin must be normal to the local surface element, it experiences a Berry phase factor originating from the 
overlap of its spin states at different space-time points. If the curvature of the surface is nonzero, such a Berry phase 
cannot be eliminated by a gauge transformation. The effect is analogous to that of a fluctuating magnetic field. This 
is precisely the origin of the fluctuating gauge field in the effective theory, Eqs. (44) and (|45|), of the Hubbard (or t-J) 
model. 

The layout of the rest of the paper is as follows. In Sec. we describe the path integral formalism for the Hubbard 
model. In Sec. [II, we apply this formalism to the ID Hubbard model, show how the standard results are reproduced, 
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and develop the analogy with Polyacetalene. In Sec. IV, we apply the formalism to the 2D Hubbard model, discuss the 
connection with the spin-bag model proposed by SchriefFer, Wen, and Zhangu, examine the quasiparticle formation, 
and derive an effective action. In Sec. 0, we give a conclusion and a summary of our results. 



II. PATH INTEGRAL FORMALISM 



The Hamiltonian of the Hubbard model is, 

H = -t J2 (cLcjT + c]^c,^) + C/Xl'^^T%-i' 



(1) 



<l,]>,<7 



where c|^^ is the electron creation operator and = The summation over <i,j> only includes the nearest 

neighbor links. The partition function is given by 



Z = 



I vie 



c\ exp 



dr 



r{dr - fJ-)Cia + H 



(2) 



where /i is the chemical potential. 

The usual way to proceed from here is to factorize the Hubbard interaction via a Hubbard-Stratonovich transfor- 
mation. There are many ways to factorize the Hubbard interaction, and in principle they are all equivalent. However, 
if one does not perform the functional integral over the Hubbard-Stratonovich field exactly, in particular if one only 
integrates over the long-wavelength and low-frequency components of the auxiliary field, then different decoupling 
scheme gives different results. The latter is a well known phenomenon. For example in the mean-field theory (i.e. one 
only keeps the q = 0,uj = components of the Hubbard-Stratonovich field) it is crucial to find the "right" decoupling 
scheme, or equivalently the right order parameter. The situation can be even worse. In some cases the low-g and 
low-w components of the Hubbard-Stratonovich fields associated with different decoupling scheme captures different 
aspect of the low-energy physics. When that happens one has to include all relevant auxiliary fields in the effective 
action. For a prior example of such multi-auxiliary-field treatment, thejxader is referred to Finkelstein's derivation 
of the non-linear a model for the correlated disorder transport problem.Ej 

For the Hubbard model near half-filling, we expect that the important charge fiuctuations are near momentum 
(f = 0, while the important spin fiuctuations in 2D are around q = (7r,7r). In the following we shall use {A}^^^^ to 
denote the slow A- modes around the modulation wave vector qg and implicitly impose the condition < tt. Inspired 
by the above understanding the Hubbard interaction factorizes into the following fiuctuation channels. 



j l9l<'r 



(3) 



where tt — (tt, tt). Note that the factor 1/4 in front of the first and third terms on the right hand side of (||) is 
obtained by omitting the following two terms in the factorization: (l/4){c|C| — c.lCl\^{c^^c^ — CjCjl-g- and (l/4){c|C| -I- 
cj^cj^lff+g- {c|C| -|- C|C|}_ff_g-. In the above we have explicitly omitted the Cooper channel. The necessary Hubbard- 
Stratonovich decoupling thus includes 



i 



Y V^c('?){cTCT +CiCx}_5-+ ^Aciq)^c{-q) 

-fA^(7r {c|C| - C|cJ_ff_g- -I- — A^(7f + (7)A^(-7r - (?) 
+A+(7f + q) {c|cJff+5-+ {qCT}ff+,-A+(7r + q) + ^\A+{Tf + q)\' 



(4) 



In coordinate-space we have 



3 



+ (-!)■'■[ A+(j) cjtc.x + c.iCjT A+O')] + ^|A+(j)P 



(5) 



In the above, we have explicitly pulled out the staggering factor from A.z{j) and A_|_(j). Therefore, all the remaining 
variables Ac(j), and A+(j) are smooth fields. In 2D, each lattice site is labelled by two integers, j = {jx,jy), 

in the unit of lattice constant, and the staggering factor denotes (—1)-' — {—ly^'^^". 
Combining Eqs. to the partition function can be approximated as. 



Z 



I?[c,c,Ae,A,,A+,A+] exp 



dT 



(6) 



H' = -t (c„c,. + /i.c.) + ^5]{[A,(jr + [A,(jr + |A+(j)p} 



<i,j>,(j 



^z{j)crz + ^A+{j){ax + i(jy) + ]^A+{j){ax - iay) 



(7) 



a/3 



The fact that only smooth configurations of Ac(j), Az(j) and A+(j) are included in the path integral ensures that 
there is no overcounting. 

To emphasize the SU(2) symmetry we group Az{j) and Aj^{j) together and write them as 



Az{j)<7z + '^A+{j){(T^ +i<yy) + ^A+(j)(cr^ - i(Ty) = As{j)g+{j)<Tzg{j). 



(8) 



Here As{j) is a real scalar field representing the longitudinal component of the local moment, while g{j) is an SU(2) 
matrix encoding the transverse components of the local moment. We emphasize that both As{j) and g{j) are smooth 
functions of space and time. 

Let us now perform a space-time-dependent local gauge transformation to absorb the g and g+ in Eq.(^). 

C^ = C^g+{i), with C^ = (C,|,C4). 

We then rewrite everything in terms of the new variables. A straightforward manipulation shows thatH 



(9) 
(10) 



Z = J V[c,c,A„A,,g] expl-^^drrj , 

^ - E + 9ij)drg+{j) + (-!)■' A,(j>, + tA,{j) - f,] c, 
j 

-t J2 [^»5(05+(j)c, + h.c] + i ^ [Alij) + AlU)] . (11) 

<i,j> j 

In the above we have omitted the primes on the Grassmann variables. Our subsequent treatments of the path 
integral ( [ll| ) involve a saddle point analysis followed by expansions that include small fluctuations around the latter. 
In specific our saddle analysis assumes time-independent (but spatially dependent) fields given by 

A,(r,j) = -*A(")(j-), 
As{r,j)^Ai'\j), 
9{tJ) = ffo(j), 
and the corresponding saddle point equations are 



(12) 

(13) 
(14) 



(15) 
(16) 

In the above rij = '^la'^j<^- This analysis is equivalent to the Hartree-Fock solution. It is also equivalent to the 
variational solution constructed using a Slater determinant as the trial wave function. 
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III. THE ONE-DIMENSIONAL HUBBARD MODEL REVISITED 



The purpose of this section is to demonstrate that our approach reproduces the well known results in ID. Hopefully 
it will strengthen the readers' confidence (as well as ours) when it is applied to 2D. 



A. ID Hubbard model at half-filling: the non-linear a model 

To warm up, let us first consider the ID Hubbard model at half-filling (i.e. with on average one electron per site). 
It is known that for arbitrary repulsion, the half-filled Hubbard model in ID is a charge insulator, and the only 
low-energy excitations are spin in nature. Here we derive the effective theory for such spin excitations. 

We first find the saddle point solution in the form of Eqs.(p^)-(p^). The result is 

A,(j)^-zA(«), A,(j) = Ai°), g{j) = l. (17) 
The corresponding saddle point (or mean field) Hamiltonian is 

(cjc, + h.c) +Y,c\ [(-l)'A(")a. -p\c. + ^ [(A(°))^ - (A^"))^] N, (18) 



Hmf — — i 

<j,i> 

where N is the number of lattice sites and 



/i = ^-A(°). (19) 



This Hamiltonian can be easily diagonalized. At half-filling /i = and aI"'' satisfies^ 

2 /"^^ dk 1 



^ ^(2tcos/fc)2+ fA,*°^ 



(20) 



For U/t = 8/3, 4, and 8, the solutions to Eq.(^) give « 0.726^, 1.539t, and 3.754^ respectively. The solution 
describes a long-range ordered antiferromagnets with the local moment ttiq ~ 0.272, 0.384 and 0.469 respectively. 
The mean-field electronic spectrum has a gap separating an empty conduction band from a full valence one. 

To study the low-energy excitations, we include fiuctuations in Ac(j), As(j) and g{i) around their saddle point 
values, and integrate out c and c. Since the action is quadratic in the fermion variables after Hubbard-Stratonovich 
factorization, the integration produces a fermion determinant. It turns out that both Ac(j) and ^s{j) fluctuations 
are gapful, thus can be safely ignored at low energies. We expand the fermion determinant to quadratic order in 
gd^g^ to obtain (the details are supplied in appendix B) 

£eff = Y.^St)T^ [a.g{i)d.g+{i)\ + ^ j dx ^ [Tr {d^gd^g+) + Tr {a.gd^g+ a.gd^g+)] . (21) 

In t he ab ove, (S'f) — {c\i7zCi)/2 = (— l)*mo is the mean field magnetization, and ps is the spin stiffness given by 
Eq. ( p2^ ). The spin wave velocity turns out to be the same as the bare Fermi velocity. The second term in ( pT| ) is 
the well known non-linear a model written in terms of the SU(2) matrix g{i). Using n ■ a = g'^azQ, one can easily 
verify that 

{d^nf = Tr {d^gd^g+) + Tr {a,d^g+a,gd^g+) . 
The first term of (|2l| ) represents a Berry phase. To see this, we express 

(22) 



-Z2 Zi 

where z*Za- — 1. In terms of z* — (z^, Z2), n — z*az. It is straightforward to show that 
mo^(-l)^Tr {a,g(i)d,g+{i)\ = 2too ^("l)'^*(05r2(«) 



i— \ dx dxfi ■ {fi X drfi), 9 = 27rmo. (23) 
27r , 
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While our approximate calculations do yield the right form for the effective|-aption, it makes a serious mistake. It is 
well known that the 9 in the effective action is quantized to have the value ttJIj Our approximate calculation gives a 
that depends on the magnitude of the local monient. It is also known that the low-energy behavior of the non-linear 
sigma model critically depends on the value of ^.Olj In particular, the spin excitation spectrum is gapless only when 
— (2n -I- l)7r. Since toq < 1/2, we would conclude that except for the case of [/ = oo (in that case mg — + 1/2) the 
spin spectrum is always gapful. This is, of course, incorrect. Then what is wrong with our derivation? 

The answer relies on the fact that the value of 6 reveals a symmetry of the problem. The original ID Hubbard 
model is invariant under the translation by one lattice spacing. Since it is known that the ground state of the model 
preserves this symmetry, it better be true that any low energy effective theory is revealing it. This is not so in Eq. (p3) 
unless 9 = integer x tt. Indeed, as a; ^ a; -|- a, toq —mo which invert the sign of 9. Since J dx d^n ■ {ft x drfijis 
an integer, the latter is a symmetry only when 9 is an integer multiples of tt. 

Thus if we did the calculation exactly, 9 should turn out to be tt. This has been shown to be the case recently by 
Nagaosa and OshikawaE3. The point is that approximations are not allowed in-extracting quantized quantities. To 
proceed, we offer the following symmetry based arguments to fix the value of 9iS Since from translational symmetry 
9 = integer x tt, our job is only to determine which integer. For this purpose we recall that as a function of U/t the 
half-filled Hubbard model experiences no phase transition. Thus it is sufficient for us to determine 9 in the U/t oo 
limit. Since in the latter limit the leading order result becomes exact, and since 9 remains quantized as we vary U/t, 
we can deduce its value exactly. As U/t — * oo, mo — > 1/2 and — )■ tt, thus we conclude that the exact value of 9 in 
Eq. ( p3| ) is TT. 

Summarizing the above discussions, we have arrived at the following low-energy effective theory for the ID Hubbard 
model at half-filling: 

S'eff = J dr J dx dxfi ■ (n x 9^n) -I- y 

The consequence of this action is that transverse fluctuations destroy the long-range order in the saddle point solution. 
The result is a quasi-long-range antiferromagnetically correlated state in which the correlation between local Neel order 
decays as the inverse of distance. 



(24) 



B. The charge and spin excitations away from half- filling 



In this section, we shall show that away from half-filling the elementary excitations of the ID Hubbard model 
are holons and spinous. In addition we shall demonstrate the similarities between the holon and the soliton in 
Polyacetylene. 

Our starting point is the saddle-point solution of Eq. (|l^) . Naively one might take the half-filled solution and take 
out an electron from the valence band. This is dangerous since one overlooks the effect of doping on the saddle-point 
solution. For example in the presence of a static hole, we do not expect the saddle point solution to have spatially 
uniform spin and charge density. This is because the spin and charge densities are self-consistently determined by the 
occupation of the electronic states. Around the hole, the electron density is reduced so that the magnitude of {Sf) 
must also be reduced. Motivated by the above considerations we look for a self-consistent solution of the following 
form: 



A,(i) = -zA(«)(j), A,(j) = A(°)(j), = 
The saddle point Hamiltonian is 

H^f = -t E (c.k+/i.c.)+Ec][A(°)(j)-A^+(-lFA?^(jV. 

<i,i> J 



(25) 



(26) 



In Appendix |^, we show thatj-the continuum version of Eq. (p6| ) is very similar to the corresponding continuum 

e.El 



Hamiltonian for Polyacetylene .lII The saddle point equations are 

A(°)(z) = -(-l)'^[(ctc,T)-(4,c,,) 



(27) 
(28) 
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We solve the above equations numerically in a finite chain of 100 sites with 99 electrons under the open boundary 
condition. Armed with the experience from Polyacetylene, we choose the starting point of the self-consistent iteration 
to be 

A(«)(*) = A(«) tanh , ^ ^, (29) 

where xq and are two arbitrary parameters and aI'^'' is the solution at half filling. It turns out that for small U/t, 
the self-consistent iteration preserves xq . The resulting self-consistent saddle point solution suggests that xq does not 
necessarily center at a high symmetry position (such as a lattice site or the mid-point of a lattice bond). To ensure 

that the 100-site chain is sufficiently long ,11 we have checked that the difference between the values of A^^'' under 
periodic and open boundary conditions at half filling only shows up at the fourth significant figures. For U/t — 8/3 
the saddle point solution and the associated spin and charge profiles with the kink centered at the middle of a lattice 
link are shown in Fig. || and Fig. || respectively. l— . 

We notice that the boundary effects die off a few lattice spacings away from the two ends of the chain.Ea The 
ground state energy at the saddle point is — 0.196i (— 0.269Ai°^) for U/t = 8/3, where the reference energy is chosen 
to be the energy of the "doped semiconductor state" obtained by removing one electron from the top of the valence 
band of the half-filled self-consistent band structure under open boundary condition. For other values of U/t ratios, 

the energies of the saddle point solutions with a;o situated at the mid-point of a bond are — 0.242t (— 0.260Ai'''') for 

U/t = 3, -0.366t (-G.238Ai"') for U/t = 4, and -0.643t (-0.17lAi°^) for U/t ^ 8. 

The one-electron spectrum for U/t = 8/3 is shown in Fig. ^. As in Polyacetylene, there are two midgap states. 
These states have opposite spin quantum numbers, and are almost degenerate. Their wave functions are shown in 
Fig. ^. All energy levels below these two states are occupied. The spin and charge densities of the saddle point 
solution can be fitted to simple analytic functions. The least-square fitting gives, for U/t = 8/3, 

i-lT{Sn = ^ tanh , 6 = 2.978, (30) 

K:)=l ^-T. = 3.209. (31) 

2^,cosh2 (i^j 

The quality of the fitting is demonstrated in Fig. ^. In principle, we still have to show that we have allowed sufficient 
variational degrees of freedom in our mean-field ansatz Eq. (|2^) . Here we skip Jtiiis issue by referring the reader to the 
similar situations in Polyacetylene where the the soliton stability was proven.ES 

Intuitively one expects that the kink energy reaches optimum if the center of the kink is located at high symmetry 
positions of the chain. It turns out that for e.g. U/t = 8/3 the kink energy has no measurable dependence on its 
position! This is, of course, what one expects based on the continuum Hamiltonian we derive in Appendix ^ The 
validity of the continuum approximation is that the self-consistent spin and charge densities have smooth spatial 
variations. This requirement will eventually be violated when U/t 3> 1, in which case the kink is a very localized 
object. In the latter case the center of the kink becomes a crucial parameter, not only does it affect the energetics 
but it also determines whether self-consistency can be reached at all. Indeed, we have found that for large U/t the 
only possible self-consistent kink solution is a kink centered at the mid-point of a bond. 

The static kink solution discussed above represents a snap shot of the charged quasiparticle of the ID Hubbard 
model. If the time scale associated with the motion of the quasiparticle is long compared with the inverse excitation 
gap of the soliton band structure, to a good approximation we can regard it as a rigid boost of the static soliton 
discussed above. In that limit, in order to determine the quantum numbers of a moving soliton it is sufficient to 
study the corresponding quantities assuming the soliton is at rest. In the following we shall assume such an adiabatic 
picture. Let us first determine the charge of a kink. Since we are interested in the case where a continuum description 
is possible, we shall concentrate on the case of U/t = 8/3. Moreover we shall assume that the center of the soliton 
is situated at the center of a bond. In Fig. ^ we show the change in site occupation number due to the presence of a 
kink. From Eq. (|3l|) we deduce that 



J2{Sn^) ^ - / dx — ^ = -1. (32) 

i J-^ 2ecCosh2(^j 

Therefore, the kink carries charge e. 

Now we determine the total associated with a kink. This is sufficient for the purpose of later discussions because 
it turns out that plays the role of coupling constant in the residual gauge coupling between the soliton current 
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and the fluctuating antiferromagnetic background (see later). The total (S^) associated with a kink is given by, 
(S^) = E,(^pair(0>, where (S^^U^)) = (^'(*) + S^i + l))/2. In Fig. § we show the {3^,^^)) Profile, where we 
see that the only nonzero (»S'pair(*)) appears near the kink. (The zigzag pattern near the chain ends is caused by the 
boundary and it is not shown in the figure.) Had we plotted the same quantity for the half-filled case we would have 
found zero everywhere except near the ends of the chain. From Fig. ^, It is clear that by summing up all {Sp^-^j.{i)) 
we get zero. 

Thus on a scale larger than the size (determined by the larger one between and ) of the kink we can regard 
it as a point object carrying charge e and zero . These quantum numbers are consistent with those of the holon in 
the ID Hubbard model. 

Next we turn to the neutral spin excitations. In the presence of the longitudinal spin order one naturally expect 
the low-lying spin excitations to be described by the smooth space-time fluctuating g{T,i). However, is that all? For 

instance can we have a neutral spin excitations with a sharp kink form in A^''^ (i) ? In the following we show that 
the answer is no. Speciflcally, we shall start from the solution corresponding to a charged soliton and show that if 
the missing electron is put back into the system the soliton decays into a smooth spin twist. For this purpose, it 
suffices to find a path connecting the configurations corresponding to the kink and the smooth twist, and show that 
the energy continuously decreases along the path. Considering Eqs. (0) and (H), the path we found is. 



+ (-l)''A(") 



1- A 



2^^cosh2 (i^) 



tanh 



I - xo 



cosh 



( i-XQ \ 



(33) 



The first segment of the path is given by increasing A from zero to one while keeping and fixed. We note that 
when A = 1, the spin order profile given by Eq. ( ^3| ) becomes a twist of the size ^g. The second part of the path is to 
increase while fixing A = 1. Again we choose a chain of 100 sites with U/t — 8/3 and xq situated at the center of a 
bond for the demonstration. To obtain the electron wave functions along the path we modify the Hamiltonian ([2^), 
by replacing the onsite potential term with the right hand side of Eq. (|33|). We diagonalize this Hamiltonian to 
construct the wave function by filling 100 electrons in the lowest 100 energy levels. Then we obtain the variational 
energy by calculating the expectation value of the full Hubbard Hamiltonian. The result is shown in Fig. ^, where the 
the reference energy is chosen to be the energy expectation of the half-filled mean-field state. In the left panel, the 
energy is shown as the function of A while fixing = 2.978 and = 3.209. We see that the kink is unstable once an 
additional electron is introduced. This is not difficult to understand if we recall the physics of Jahn- Teller effect. The 
two nearly degenerate midgap states are empty in the case of charged soliton. If the missing electron is put back, it is 
energetically favorable to split the two midgap levels and occupy the lower level. The off-diagonal part in Eq. (^3) 
precisely provides matrix element between these two states and consequently split them. In the right panel of Fig. 
the variational energy is shown as a function of the twist size while fixing A = 1. The twist given by Eq.r(g3[) with 
A = 1 further decays to an infinitely smooth one, as illustrated in Fig. ||. In fact, the twist is just a spinoni 



C. Effective theory 



In the above discussions the saddle point solution for doped ID Hubbard model consists of static and localized 
solitons (each accommodating one hole) separating the otherwise perfectly ordered antiferromagnetic domains. Ig- 
nored in this mean-field theory are the smooth space-time variation of the direction of the order parameter, and the 
propagation the soliton configurations. In general it is not clear that the saddle point solution in the presence of N 
holes is N separate solitons (for example, it is important to compare the mean-field energies associated with, say, a 
phase-separated solution with that of N equally spaced solitons). In the rest of this section we will assume that the 
individual soliton remains stable against finite doping. Under that assumption we will show that the low-energy theory 
is precisely the effective theory for the Hubbard model obtained via non-Abelian bosonization (see Appendix 

We now briefly describe the derivation of the effective action, and leave the details in Appendix ^ The basic 
strategy is to allow the solitons to have time-dependent locations, and twist the spin directions in a space-time 
dependent way. Then we evaluate the cost in action due to such distortions. 

At the length scale longer than the kink size ^s, the kink (^^ can be simply viewed as a point at which the 
longitudinal spin order parameter As{i) flips its sign. Consequently at a given site the sign of As{i) has changed as 
many times as the number of holons to its left: 
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A,(a;) ~ A'"' cos 



where 



is the holon density, 
small perturbations: 



n{y) 



TT / dyn{y) 



5{v-Xi), 



0<x<L, 



E 



(34) 



(35) 



At low energy and long wavelength we expect the holon density to be its average value plus 
n{x) — Ufi + dn{x). The spin order profile in the finite doping becomes, upon replacing the 



lattice site j with the continuum coordinate x, 

(-1)^A,(.?) 



"'0 



A^'') cos [2kFX + ^ix)] 
dy Sn{y). 



where kp = 7r(l — n/,,)/2. Furthermore, the phase ^{x) satisfies the following density-phase conjugation relation: 

d<^{x) 



dx 



— — TT 6n{x). 



(36) 
(37) 

(38) 



Substituting Eq. ( |36| ) into Eq. (|l 1\) , and integrating out c and c, we derive the desired low-energy effective theory in 
the case of finite doping. The charge part is given by a free boson theory 

r0 /• r 1 



c(c) _ Pc 



dr / dx 



(39) 



wherepc and Vc are the charge stiffness and velocity, given by Eqs. ( ^29|) and (|B30| ). The spin part is still given by 
Eq. (p^). Both Pc and Vc are different from ps and vp. The detailed calculation is included in Appendix]^, (also see 
Reference^.) |-| 

As noted by Lee from the Bethe Ansatz solutionQ, the Fermi velocity in the low-energy effective action is not much 
different from the bare one of the non-interacting case. This is consistent with the fact that in Appendix ^ the 
integral in Eq. (B26) is independent of A^^\ We note that the parameters in the effective action are subject to the 
renormalizations due to short distance fluctuations which are simply dropped when we take the continuum limit in 
Appendix However, as the order of magnitude is concerned, the results are reliable. This is in sharp contrast to 
the situation in two dimensions where the effective Fermi velocity is renormalized down to the order of J^r^ /U), as 
we shall see below. 



IV. THE HUBBARD MODEL IN TWO DIMENSIONS 

Since the discovery of the high temperature superconductors, AndeisonB has been suggesting that the one-band 
Hubbard model in two dimensions (or the closely related t-J modelO) captures the essential physics of the Cu02 
plane in the cuprate superconductors. Ever since then, the properties of 2D Hubbard model (on square lattice in 
particular) has been a focus of investigations. An important parameter for the cuprates is the doping level. For 
example, within a class of compounds, the samples with the highest Tc usually have doping level around 15%. In 
terms of the Hubbard model it means that the averaged site occupation number is around 0.85. As in one dimension 
the starting point of our discussion is the half-filled Hubbard model. 

By now it is a common consensus that at half-filling, the Hubbard model describes-an ordered antiferromagnetic 
Mott insulator and the long wavelength effective theory is a 2-1-1 D non-linear a modeled. The derivation of the latter 
from the Hubbard model proceeds in the same way as in one dimension except that the Berry phase term is canceled 
out in two dimensions O The antiferromagnetic order present in the half- filled 2D Hubbard model reflects the fact 
that the non- linear a model can be ordered in 2-f 1 dimensions. 

An important question is how to proceed when the antiferromagnet is doped. In the conventional approach, one 
performs the Schricffer-Wolf transformation on the Hubbard model to eliminate states having doubly occupied sites in 
the U /t oo limit. The result is the familiar t~J model. In the following we shall adopt the same approach as the one 
used in the last section. As we shall see, there are two levels of quasiparticle formation. At the first level, like in ID, 
we shall find that in the mean-field theory doping in the antiferromagnetic long-range ordered state produces solitons, 
whose corresponding band structure exhibits states inside the gap. Unlike ID, the solitons are non-topological. At the 
second level, due to the residual coupling between the soliton and the spin degrees of freedom, the soliton is further 
dressed so that the bandwidth associated with its hopping is reduced from ~ t to ~ J. In ID the second level of 
renormalization is absent. 
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A. Mean-field solution for one hole 



The Hubbard-Stratonovich decoupling scheme and the subsequent mean-field ansatz are the same as in ID. In 
particular, the mean-field Hamiltonian is 

Hmf^~t ^ (4c,+/i.c.)+5]c][A(«)(j) + (-l)^-^+^='A(°)0V,]c, 

<i,j> j 

+ i5:j[A(")0-)]'-[A(°)(,)]'|. (40) 



The only difference is that now i,j label the square lattice sites. The corresponding saddle point equations are similar 
to Eq. @ and 



U 



(41) 
(42) 



We solve the above equations numerically on a 14 x 14 lattice with periodic boundary condition. The number of 
electrons is such that there is only one hole. We have chosen three different starting points for the self-consistent 
iteration. In each case, we reduce the Neel order parameter from its value at half-filling symmetrically around a lattice 
site, lattice link, and the center of a plaquette respectively. For U/t not too small (say U/t > 10/3), the first two 
starting points converge to a non-topological soliton centered on a lattice site. We have also checked that a small 
perturbation in the starting configuration does not affect the final self-consistent solution. The magnitude of the 
resulting order parameter has a fourfold rotational symmetry about the center. For U/t > 10/3, the third starting 
point first converges to a diagonal cigar-shaped non-topological soliton situated at the plaquette center. However, 
if we slightly break the refiection symmetry about the plaquette center along the ridge of the diagonal cigar-shaped 
spin profile, the iteration further converges to the same solution as found from the previous two starting points. For 
U/t > 10/3, the self-consistent mean-field solutions are non-topological solitons. The typical spin and charge density 
profiles are shown in Fig. and the corresponding mean-field band structure is shown in Fig. We note that two 
midgap states are present. For one hole, only states in the lower band are occupied. The energies of the soliton with 
respect to a hole at the top of the valence band of the half-filled band structure are SE = — 0.051t (— 0.049Ai''-') for 
U/t = 10/3, -0.386^ (-0.155Ai°^) for U/t = 6, -0.650t (-0.182Ai°^) for U/t = 8, -0.858t (-O.ISSAI"^) for U/t = 10, 
and — 1.017< (— O.lTQAi*^-*) for U/t = 12. For U/t > 8, the calculations arc done in 10 x 10 lattice since the finite size 
effect is small for large U/t ratios. 

In Fig. |l^, we show the changes in the spin density induced by a soliton. We have checked numerically that 
Ei a ^(cLc^a) = -1 and I Y,i \ = 1/2. Thus, unlike ID, the soliton carries both charge and spin. 

The same mean-field studies for one hole in the 2D Hubbard model has been carried out by Su and Chen,El and 
Choi and Mele.El For large U/t our results are consistent with the their findings. For example, we have found that 
for U/t> 10/3, the saddle point solution is a soliton whose associated spatial variation of the Neel order parameter 
respects the fourfold rotational symmetry of the underlying lattice. Moreover, our self-consistent band structure 
agrees with that found in Ref. ^ (for U/t = 5 on 10 x 10 lattice) within numerical uncertainties. However, for smaller 
U/t (say U/t = 2) we found that unlike previous claims that the solitons are cigar-shaped, after a large number of 
self-consistent iterations the soliton converges to a linear bag running diagonally across the whole finite lattice. Based 
on the length scale defined by z)f/2A1°^ we expect significant finite size effects for small U/t. In addition, for small 
U/t the detailed structure of the non-interacting band (such as the magnitude of the next nearest neighbor hopping) 
starts to have a stronger effect on the soliton shape. At present we have not studied enough lattice sizes and different 
boundary condition to establish that the cxteiadcd soliton is the genuine solution in the thermodynamic limit. In that 
context it is interesting to observe that SchulzC3 suggested that at finite hole concentration, the carriers segregate into 
linear walls separating 7r-phase-shifted antiferromagnetic domains. We take the appearance of one-hole line-shaped 
soliton for small U/t in our study as an indication in favor of the former possibility. 

The saddle point solution for large U/t is robust and simple to understand. Numerically we have found that when 
U/t is increased, the size of the soliton is reduced. Since the soliton has the same quantum numbers as a bare hole 
in the t-J model, we interpret it as a finite U analog of the latter, which, among other things, also implies that the 
mean-field solitons obey Fermi statistics. In the following we shall show that the low-energy fluctuations around the 
saddle point solution is precisely the low-energy dynamics described by the t-J model: 
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<i,j> <id> 



(43) 



In the context of our above discussions the perfect Neel state is the analog of our half-filled mean-field solution, and 
CiCT acting on the latter produces our soliton. The remaining dynamics (including the hopping of the soliton and the 
flipping of the spins) described by Ht-j are captured by fluctuations around the saddle point. 

The derivation of the eS'ective action is similar to the ID case. However, unlike in ID, the 2D solitons do have 
preferential positions (i.e. the sites) even for small U/t. Thus instead of a continuum description, we use a tight- 
binding language to describe the soliton hopping. The resulting effective action contains a 0(3) non- linear sigma 
model part describing the antiferromagnetic fluctuation, and a tight-binding spinless fermion part representing the 
soliton hopping on the lattice. Finally, unlike in ID, the soliton motion docs couple to the a model fluctuation. 
This effective action is precisely the coherent-state functional- integral action of the t-J model in the slave- fermion + 
Schwinger boson representation. The details are supplied in Appendix M. 



B. Dressing the non-topological solitons 

In the above we have argued that the effective Hamiltonian governing the motion of the solitons and the antiferro- 
magnetic fluctuations is the t~J model. Due to the residual coupling between spin and charge, the soliton motion in 
2D is a highly nontrivial problem. A lot of work has been done on this subject in this context. In particular, two of 
us have recently presented an intuitive picture that explains how the frustration of hopping in an antiferromagnetic 
background renormalizes the hole bandwidth from ^ t to ^ jJHiJ 

The question of charge motion in the presence of finite doping concentration is a much more subtle issue. This is 
due to the following facts, i) At finite doping it is not clear what is the best saddle-point solution. For example is 
the anti-phase domain wall solution suggested by Schultz the correct mean- field solution for certain range of U/t ? 
ii) Even one assumes that the individual bag-like solitons remain the lowest-energy mean-field solution, one still has 
to face the fact that in 2D, due to the residual coupling between the solitons and the antiferromagnetic fluctuations, 
the solitons are further dressed into quasiparticles. When the size of the one-hole quasiparticle reaches the inter-hole 
spacing the meaning of a quasiparticle is no longer clear. Moreover, the presence of finite carrier concentration may 
feed back to qualitatively change the low-energy dynamics of the antiferromagnetic fluctuations. Thus even though 
in the case of one hole, one has reasonable confidence that the bag-like soliton is first formed and subsequently is 
dressed into aj-quasiparticle, the same conclusion should be subject to critical scrutiny in the case of finite doping 
concentration.^ 

With these caveats in mind let us imagine that the doping level is relatively low so that we can view the finite 
doping as a collection of one-hole quasiparticles. Let us imagine that the magnetic long-range order is destroyed. 
Let us use|-tiie antiferromagnetic fluctuations with wavelength shorter than the magnetic correlation length to dress 
the solitonnJ. What remains is to determine the effective theory governing the motion of these dressed objects that 
includes their interaction with the remaining long wavelength magnetic fiuctuations. 

Since the residual long wavelength spin excitations correspond to smooth fluctuations of the spin directions we 
describe it by a non- linear a model (see Appendix |^ for derivation) , 

Snism = y / drld^ Tr {d,gd^g+ + a,gd^g+a,gd^g+) , (44) 

where ps is the stiffness and is the spin wave velocity. In particular, we choose ps so that the resulting a model 
is in the disordered phase. (Here we should point out the fact that how to get a spin liquid described by the 
disordered a model from a microscopic spin model on a lattice is still an open issue.) To find the effective action 
for the quasiparticle propagation, we need to generalize the quasiparticle construction in Ref. ^ to allow a smoothly 
fluctuating spin background. Since the dressing envisioned in Ref. ^ is a short wavelength and high-energy process, 
it will not be significantly affected as long as the spin correlation length and correlation time are sufficiently long. 
The generalized construction of the quasiparticle for the t-J model is described in Appendix |b| where the soliton 
hopping is treated exactly while the spin-exchange part is treated on average. The quantum spin fluctuation generates 
hopping of the dressed object as a whole. Let i and j be two lattice sites with quasiparticle creation operators f} and 
/J. The hopping matrix element between these two states is due to the spin exchange part of the t-J Hamiltonian, 

Hj — J^^ij^ Si ■ Sj. With only short range order, the local spin directions at the sites i and j are not necessarily 
parallel to each other. Since the quasiparticle spins are constrained to follow the local spin directions, the two 
quasiparticle states at the sites i and j do not have parallel spin directions. As a result, the hopping amplitude 
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acquires a Berry phase which is expressed in terms of the SU(2) matrix g{i). The details of the derivation are 
included in Appendix Taking into account this additional factor, the quasiparticle hopping is described by, 

ieA.r=V2,2 ieB,r=V2,2 

+ E A(r){[5(*)5+(z + r)] ^2 + /i.e.}, (45) 

where the matrix subscript denotes the respective element, and the notation r — \/2, for example, indicates summing 
over all neighboring sites of V2 units of lattice constant. The first two terms of Eq. (|4^) represent the hopping among 
the same sublattice sites, while the last one describes the quasiparticle hopping term between different sublattices. 
Long range hopping farther than three lattice spacings has been neglected. We expect that the results for the 
coefficients a(r) and X(r) calculated in Ref. ^ remain a good estimate so long as the correlation length (or time) of 
the antiferromagnetic order remains long. The two parameters, a{V2) and a{2), describe a quasiparticle ha n d Q fj-a. 
width of order 2 J, with the minimum at (7r/2,7r/2) and extended van Hove regions around (7r,0) and (0, tt) j3~113'ElrE£l 
Alternatively, these parameters, as well as A(r), can also be taken as phenomenological parameters of the effective 
theory to be determined from fitting experimental data. 

Using the representation (^) of the SU(2) matrix, Eq. (^) can be recast as 

= E "WE^^(^)^<^(*+^)-^i^-^*+''+ E "WE^'^'^^)^^(^+^)-^i^'^*+'- 

ieA,2,2 a i£B,r=^/2,2 

+ E Xir){[zm4U + r)-z;it)zl{t + r)]fjf,+r + h.c.}. (46) 

i€A,r=V5, 3 

The corresponding part of the action is 



^qp= / '^^l l^hdrh + H,p\. (47) 
The complete effective action is 

S'cff = Snlsm + Sqp. (48) 



C. Closing remarks 

What can be said about the nature of the carriers after being dressed by both short and long ranged magnetic 
fluctuations? For one hole in an antiferromagnet it is widely believed that the soliton is dressed 'ui\£s_a, quasiparticle. 
Here we quote some numerical results. For the t-J model there exists direct diagonalization results.Bilj The conclusion 
of these studies paints a picture in favor of a quasiparticle with internal structure. These results can be understood in 
the string picture developed in Ref. For the 2D Hubbard model, most numerical results are obtained from Monte- 
Carlo simulations .ElrE3 Unlike the direct diagonalization, these results are for non-zero temperatures. Moreover, it 
is often hard to obtain very low temperature results. In any case the state-of-the-art low temperature Monte-Carlo 
results awpq sfttith the t-J model diagonalization. The outcome is a quasiparticle obeying a dispersion relation best 

fltted tMWm 

tj: = Aa{y/2) cos kx cos ky + 2q!(2)(cos 2kx + cos 2ky). (49) 

The two parameters a{\/2) and a{2) have also been calculated in the string picture.0 This result compares very 
favorably, jssdth the Angle-Resolved Photoemission results of the high temperature superconductors near optimum 



lavoraaiy. sat 
doping.ll3 a 



In the presence of a finite concentration of holes, it is far less clear whether the quasiparticle picture is correct. The 
effective Hamiltonian, given by Eqs. (^) and (^Sj), describes the ultimate dressing of the quasiparticles which have 
already been partially dressed by short wavelength antiferromagnetic fluctuations. The best chance for the relevance of 
the Hamiltonian is when the magnetic correlation length (time) is long compared with the inter-hole spacing (iniierse 
intra-quasiparticle excitation gap). This Hamiltonian has been postulated and studied by a number of authors .EZTEEI 
Here we have presented a derivation of it. 
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V. CONCLUSION 



In this paper we have demonstrated that the two low-energy dynamical degrees of freedom in doped Hubbard model 
are a) soliton charge carriers, and b) smooth antiferromagnetic fluctuations. We have arrived at these excitations 
starting from the SDW mean-fleld theory. The quantum number of the charge soliton distinguishes one dimension 
from the rest. In ID the charge solitons are the anti-phase domain wall of the longitudinal spin order. They carry 
no (5*'^), consequently they do not suffer from further dressing by the magnetic fluctuations. Thus in ID the final 
effective theory is a decoupled charge and spin model. In two dimensions, the solitons are bag-like objects and they 
do carry (S^). Consequently, they are further dressed by the antiferromagnetic fluctuations. The ultimate effective 
theory consists of two types of charged particles, i.e. with opposite fictitious charge, coupled to fluctuating gauge 
fields. The final zero temperature state of such model is currently not known. 

At last, it is important to point out that in two space dimensions we made an important assumption at finite doping 
— that individual solitons are still the lowest-energy mean-field solution. If that turns out wrong, and if the charged 
object has the form of an extended domain wall, then the discussions presented here are irrelevant. 
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APPENDIX A: TAKAYAMA-LIN-LIU-MAKI EQUATION 



The-one-dimensional Hubbard model written in the form of Eq. (|11|) does not resemble the Su-Schrieffer-Heeger 
modeim describing Polyacetylene in the lattice even if the transverse fluctuations of the spin directions are turned 
off. One difference is that the longitudinal spin order As{i) does not have dynamics in Eq. (pT[). In other words, 
Eq. ( pi] ) does not contain time derivative of As{i). This difference is not so important since the quantum fluctuations 
will generate dynamics for As(i). Another difference lies in the way how the longitudinal spin order is coupled to the 
electrons. Despite these differences, we show in the following that the continuum version of the saddle point equations 
corresponding to Eq. ( p6| ) are essentially the same as the Takayama-Lin-Liu-Maki equation^ determining the soliton 
solution in the continuum for Polyacetylene. 

In the continuum limit, we substitute into Eq. (Eq) 



(Al) 



where iPl,r{''') are slowly varying fields in space, V'I(j) = (^I|(j): V'Ij^(j))i kpa = tt/2 at the half filling, and a is the 
lattice constant. The saddle point Hamiltonian, Eq. (p6|), becomes 



(A2) 



+ 



Y,Af\j) [v^i(j>,v^fi(j)+^is(j>.^L(j)] +^e{[^^°^w]'- 



The discrete summation can be converted into integration according to 



(A3) 



The two spin components are decoupled in the above mean field Hamiltonian. The mean field Hamiltonian can be 
rewritten as 



mf 



E 

Q = ± 
1 



dx 



dx 



- A(°)(x)' 



(A4) 
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where vp = 2ta, and Ux are Pauli matrices. The spin index a = ± corresponds to spin up and down respectively. 
This continuiLfa version of the saddle point Hamiltonian has the similar form as the continuum Hamiltonian for 
PolyacetyleneBZl Expanding the operators ^l.r{x) in terms of their eigenfunctions 

V'LaCa;) = ^ u„^Q(a;)V'LQ(K), (A5) 

tpRa (x) ^^Vi^^a {x)i^Ra (k) , ( A6) 

K 

the corresponding saddle point equations are very similar to Takayama-Lin-Liu-Maki equation]^ 

K,.a 

(x) = eK,aVK.a{x), (A8) 

^^°Hx) = - ^ 51 H " K,a(a^)^'K,aW +<,a(a^)"K.a(x)] , (A9) 
a — ± /^Goccup 



A(°H^) = ^E E [KoWP + Ko(^)l']- (Aio) 

a— ± At^^occup 

The only difference is the appearance of an additional self-consistent equation for A^^'* (x) . This analogy underlies the 
common physics of the holon in ID Hubbard model and soliton in Polyacetylene. From this analogy, we identify the 
counterpart of the dimensionless electron-phonon coupling constant Xe-ph of the Polyacetyleneli3 to be 

27rvFXe-ph " ^' ^ ^'-"""^A^f ^^^^^ 

In the absence of the self-consistent equation for Ac"'' (x) , the saddle point equations have been solved in closed form 
giving rise to Af'\x) = Ai""* tanh(x/^). The electron density calculated using the eigenfunctions in the presence of 
the kink has a shallow dip at the kink position which has the form 1 — l/[2^cosh^(x/^)]. This shallow dip causes 

additional scattering effects when there is an additional equation for Ac""* (x) . Although the true solution in this case 
will be different from that in Polyacetylene, we expect that the difference is only quantitative. In particular, the 
change should not be significant when ^ ^ 1, i.e. for small U/t. 

The presence of the Hartree term distinguishes the continuum equations in this case from those in Polyacetylene. 
However, the localized states are still degenerate in the spin index. To see this, let us define 

fi^Hx)=u^{x)±v^{x). (A12) 

The coupled equations for Ua{x) and Wq(x) (the index k is omitted) become 

- z [vpdx ± aA(0)(x)] /(±)(x) - [e„ - A^^^x)] fi^Hx). (A13) 
We see that the equations are invariant under /j^-* /j^-* and /j^' fl^^- 

APPENDIX B: CALCULATION OF THE FLUCTUATIONS AROUND THE SADDLE POINT 

The general formalism of calculating Gaussian fluctuations is as follows. Substituting into Lagrangian (|ll|) the 
expansion 

A,(j) = A(«)(j) + <5A,(j), (Bl) 

A,(j) = -zA(o)(j)+5A,(j), (B2) 

we separate the full Lagrangian, of the form of Eq. ( [TTI ) into a mean- field part and a fluctuation part. 

C = Lra! + 5/:, (B3) 

= E E + '^{^L},!.,.}) ''^-^ + ^^^)' 

(5/:-EEc-/^ W'\SA,,SAs,g)] c,,,. (B5) 



14 



The timc-indcpcndcnt matrix ft,'^"^ can be specified once the saddle point solution is obtained. The matrix h^^^ depends 
on the fluctuating fields dAc{j), SAg{j) and g. The quadratic part, 

Co{SA,,6A,) = (1/f/) ^[[SAs(i)f + [6Ac{i)f} , (B6) 

i 

is unaffected by the procedure of integrating out fermions to obtain the effective action defined through 

The integration is carried out mainly by expansion. The first step of the calculation is to diagonalize Cmf- This is 
achieved by a linear transformation, 

K 

To have a Jacobian equal to 1 or to preserve the proper commutation relation {tpl^^ V'«2,ct} = ^ki,k2 operator 
language, the transformation coefficients must satisfy 

(B9) 

n 

In terms of the new Grassmann variables V'k.o- and V'k,(t, the mean-field Lagrangian is diagonalized, 

^mf = ^i)^A9r + e^,a)^n,a + Co{5A^,dA,), (BIO) 

where £^,0- is the eigenvalue of the matrix h^^\ Then, we can represent the interaction part of the Lagrangian in terms 
of the new variables V'k.o- and V're,o-, 

S^=J2 E^'^-mW hi^lir) ^.,,.(r), (Bll) 

Kl,K2 11,^ 

where 

hi'^:l(r) =Y.ai^> \h^'\6A,,5K,g)\ (B12) 

In the final step, we expand the interaction 5C in the partition function. To the first order, 

{5C) = ^ Y.^t'':X{i'>^.Ar) V^.„.(r)). (B13) 

In the second order of the expansion in 5C^ the generated effective action is 

^iff = E Y.^,A<^U ^2,Vn) hi%\iUr,) h^:;':l{-Un), (B14) 

where Vn = 'i.im/v is the Bosonic Matsubara frequency, and h\tuXiyn) is the Fourier transform of /(.^'^'.^^(t). The 
function n^_j,(Ki, K2, 1'n) has the usual Lindhard form 

f>.,,.{Kl,K2,iyn) = . ^ —, (B15) 

where /(e) = 1/(1 + e^'^) is the Fermi-Dirac distribution function. Since the particle-hole excitation in the spectrum 
of the mean-field state has an energy gap, we can safely approximate II^^^{ki, K2, Vn) — n^,i/(Ki, K2, 0). In the rest of 
this section, we shall apply this procedure of obtaining effective action to both ID and 2D Hubbard model. 
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1. One-dimensional case 



Substituting Eqs. @ and (|b|) into Eq. the full ID La grangian is written as 
C^Y.'''j[^--^^' '5A,(j) + + A(°) cos {2kFRj + fx. 



A 



(0) 



<ij> 



(B16) 



In one dimension, we can simplify the Lagrangian by taking the continuum limit. The fermion fields can be linearized 
around the Fermi surface, which consists of two points in one dimension. This is done by Eq. (Al) but with kp = 
7r(l — nh)/2 in the general case. Since and g{i) are slowly varying functions, we can expand 



g{i)g^{j) - 1 = ag{i)d,g+{i) + -g{i)dtg+{t) + • • • , 
where a is the lattice constant. Using this expansion, we find 

-t J2 [c^ 9ii)9^ij) cj + h.c] = -VP Y {^Aij) [id. + g{j)td.g+{j)] Mj) -{l^R}] 

<id> j 
j 

where we have denoted 

vp = 2tasm{kpa). 

This generalizes the definition of vp to the arbitrary electron filling factor. Similarly, we have 

aI"^ cos $j \i^l{j)(Jzi^B.{j) + h.c] , for kpa = 7r/2, 



A(°' Y ^3 cos {2kpRj + $j) cr^ 



e'"^'^lU)<^z'fpRU) + h.c. , foi kpa ^ n/2. 



We also note that 



2 Y cos\2kpR., + ^,)=N + Y cosiAkpR, + 2$,) ^ | ;;^V E.ts(^$lf ' 



for kpa — Tr/2. 



(B17) 



(B18) 



(B19) 



(B20) 



The other terms in Eq. (B16) are simplified in the same way. In the case of finite doping, kpa ^ 7r/2, we can redefine 
the fermion fields to absorb the phase e'*^ accompanying aI^-* in Eq. ( B19D , 



(B21) 



Substituting Eqs. (B17) through (B21) into Eq. (B16), we obtain for the case of finite doping, 



-mf 



dx 



.(0) 



2U Ua 



dr~ivpdx iAi''''o-2 
iAi^'a^ dr + ivpd^ 



dx[SAcix)Y 



6C = dx 



ijl{x)AL{x)':pLix) + 'iPI{x)Ar{x)^Pr{x) 

where we have converted the discrete summation into integration according to Eq. ( [A3| ) and 

gii)^gix), SA,{i) ^ dA,{x). 



(B22) 
(B23) 

(B24) 
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We have also introduced the notations, 



(x) = g(x){dr T ivFdx)g^ix) + -{±dr - ivFd^)^{x) + i(5Ac(a;). 



(B25) 



It is now straightforward to calculate the Gaussian fluctuations according to the general prescription, Eqs. (B4) and 
(B5). We no te that the site index in Eqs. (B4) and (^5| ) corresponds to the chirality index and the space coordinate 
in Eq. (B22). The lattice summation corresponds to 



chirality— L,i? 



dx. 



After Fourier transformation to the momentum k space, Eq. (B22) is easily diagonalized by a Bogoliubov transforma- 
tion. For 5C given by Eq. (B23), {5C) = 0. In the second order contribution from 5C expansion, the transverse spin 
fluctuation part is decoupled from the rest and is given by the non- linear sigma model, Eq. (|^. The spin stiffness in 
Eq. ( p^ ) is given by 



Ps = 



dk 
2^ 



(■0[(fc, W„)V';?(fc, W„)) 



24 

/3 



dk 
2^ 



27r' 



(B26) 



Note that the Gaussian fluctuations only give rise to the second term of Eq. ([2^). The flrst term is due to the 
topological effect in Eqs. (B22) and ( |B23| ) which is not captured by the expansion in SC. When tracing out the 
fermi ons, t he res ultin g effective action contains an imaginary term known as the Berry phase. In the language of 
Eqs. (B22) and ( B23 ), the Berry phase appears in the form of chiral anomaly. Mathematically, the chiral anomaly 
exists because there is no regularization which preserves both gauge and chiral invariance. Physically, this Berry phase 
term is the electric fleld expuerienced by the phase of the holons $(2;). This Berry phase has been calculated exactly 
by Nagaosa and OshikawaS which is given by Eq. ( p3| ) with 9 ~ n, independent of doping. We shall not repeat the 
calculation here. 

The rest of the effective action including up to the second order contributions is 



dr dx ■ 



2i 

VF 



(5Ac(a;) 



(B27) 



where pi"' = ps — vf/^tt. The SAc{x) fluctuation is gapful and therefore can be integrated out. The remaining low 
energy mode described by ^(x) has an effective action 



5cff(*) dT j dx 



where 



\ Ua + ttvf ) 27r \ 



(9x<i>) 



Va 



Ua + ttvf 



Vc = VF\ 1 + 



Ua 



Ua + irvF 



(B28) 

(B29) 
(B30) 



.(0) 



The increase of the 
The same is true in 



We observe that these coefficients, as well as ps, are independent of the expectation value A 
charge stiffness by renoripaiization indicates that the interaction between the holons is repulsive 
the doped Polyacetylene.E3 

At half-flUing, both Eqs. ( B19| ) and (B20) give different results than at flnite doping. We see tha.t a pinning 
potential proportional to cos 2$ has appeared, leading to an excitation gap for the bosonic fleld $(a;).E3 Thus, only 
spin excitations are gapless at half-filling. We recover the non-linear sigma model as the effective theory of the ID 
Hubbard model at half-filling. 

Let us return to the general case of flnite doping in one dimension. We have to further show that the various 
correlation functions can be correctly reproduced using our effective action. Since the effective action we have derived 
has a form equivalent to the non-Abelian bosonized action, we only need to show that the physical fields entering our 
effective action have the same meaning as in the non-Abelian bosonization. For this purpose, we can couple external 
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fields to the spin and charge densities in the original ID Hubbard model. The rest is to show that the generated 
external field coupling terms in our effective action have the same forms as those in the non-Abelian bosonization 
action. 

Physically, we have shown in Eq. ( ^8|) that dx'^ has the meaning of holon density. Since the location of the holon 
coincides with the hole, 9^$ has the meaning of charge density. The size of the holon is irrelevant to the long 
wavelength physics. Formally, we can also prove that dx'^ has the meaning of charge density. We introduce an 
external potential hc{x,T) and couple it to the electron density in the original Hubbard model, ^ hc{i,T)cl^Cia- 



After linearization using Eq. (Al) and phase transformation in Eq. (B21), it becomes 



dxh,{x,T) U[{x)^l{x) + i^^ji{x)i;R{x)] + / dxh,{x,T) U^j^{x)^Lix)e^'^^''+'^^^^ + h 



(B31) 



Including the expression (B31) in Eq. ( B23| ), the first term of ( B31 ) adds the potential hc{x, t) to ^{±dr~ivFdx)^ix) in 
Al.r{x)- When we integrate out the fcrmionic fields in Eq. (B22), the mixing terms of type ihc{x, T)(±dr~ivFdx)^{x) 
are generated in the second order. This gives rise to an additional term 2vphc{x, T)dx^{x). Combin ed wi th the proper 
prefactor, the final term in the effective action is indeed hc{x, T)dx^{x)/TT. For the second term of (B31), the bilinear 
fermion operators can be replaced by the saddle point average. This leads to 2{^|;]^{x)^|JL{x)) J dxhc{x, r) cos[2A:^-x + 
4>(a;)]. This term will give rise to the oscillating part of the density correlation function. 

To show that the vector field n in the effective action represents the spin direction, we couple a magnetic field h{x, t) 
to the spin operator in the original Hubbard model in the form of h{i,T) ■ cl^(Jaf3Cii3/2. Introducing this couphng 
term into Eq. (^, the last term of Eq. (Q) is changed to J^i afi '^la [Qap{i) + ha/j] Cifj + /i.e., where ha(i = {l/'i)h ■ Uafj 
and Qapii) = (— l)'As(i) z9{i)\aSi ■ We can restore the original form for the last term of Eq. (j^) by a shift of 

the integration variable: Qaji Qapii) + The modification is now transferred to the second term in Eq. (|^) 

which now becomes (l/U) TrjQ + h\'^. The mixing term in the expansion of the square of Q + h is the only relevant 

magnetic field coupling term which has the form (1/?/) Tr[Q(«)/i(j) + h.c] = (l/?7) l)*As(i)rii • hi if we use 

the fact that g~^{i)azgii) = Hi ■ a. Upon substituting (— l)*As(i) = Ai°^ cos(2A;F-Ri + into the above relation, 
we obtain (Ai^^f/) J2i cos(2fcF-Ri + ^i)ni ■ hi. The expressions for the charge ancLspin density operators and their 
scaling dimensions should be compared with the non-Abelian bosonization results.CJ 



2. Two-dimensional case 



In two dimensions, the mean-field Lagrangian, the transformation coefficients a^'^j and the eigenenergies e^.o- are 
only obtained numerically. In practical calculations, they are all chosen to be real numbers. The saddle point 
Hamiltonian is given by Eq. (^). The corresponding matrix elements are 



i<°) 



-dix,±y 



(B32) 



where i,j label the lattice sites, and /i, i' are the spin indices. The typical Ai^\i) and A'^^i) for the 14 x 14 lattice 
are shown in Fig. |l^ which is the solution of Eqs. and (|4^). The transformation coefficients a'^^j are just the 



eigenvectors obtained by numerically diagonalizing h 



Lagrangian, Eq. (B5), is identified to be 



(0) 

{i,At},{j,i^}" 



The matrix element in the interaction part of the 



h^'\5A,{j),SAsU),g) =Si,, [g{j)drg+{j)]^^-t ^ Si.,,r [g{l)g+{j) - l]^ , 

+Si,j 5^^, [i 5A,{]) + sgnfi {-ly^+^y 6A,{j)] , 



(B33) 



where 6Ac{j), 5As{j) and the SU(2) matrix g{j) are the fluctuating fields. 

The first order term in the expansion of the effective action, Eq_.(^T^), does not contribute anyt hing interesting 
since we know that the Berry phase term is canceled out in 2D.c3 In the second order, Eq. (B14), we obtain an 
action in the variables 5Ac{j), SAs{j) and g{j)- Unlike in the ID calculation, the coefficients in the effective action 
have spatial variation since we are perturbing around a saddle point solution which is not uniform in space. The 
translational invariance will be restored when the soliton motion is included. Thus, we can simply perform the spatial 
average of the coefficients in S^.^ . Because the soliton has a finite spatial size, although small for U/t ^ 1, the 
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effective action S^g contains nonlocal terms, such as J2i i' /(* ~ «')Tr [g{i)drg'^{i)g{i')dTg'^{i')], where the function 
/(r) decreases exponentially with the distance r. In the momentum space, this nonlocality corresponds to a momentum 
dependence of the coupling constant which is the spatial Fourier transform of f{i — i'). If only the long wavelength 
behavior is concerned, this nonlocality can be neglected by approximating ^, f{i — i')TT [g{i)drg~^{i)g{i')drg^{i')] — 
[^-, /(i — i')] Tr [(7(i)9T-(7+(i)(7(i)9T-g+(i)]. This is equivalent to neglecting the momentum dependence of the 
coupling constant and taking its value at the momentum k = 0. After the two steps described above, i.e. averaging 

(2) 

out the spatial nonuniformity and neglecting the nonlocality of the coefficients in S^g , the g{j) fluctuations are 
decoupled from SAc{j) and SAs{j). Both 6Ac{j) and 6As{j) fluctuations are gapful. So they can be neglected. The 
g(j) fluctuations are described by the non- linear sigma model, Eq. (55), with the stiffness constant given by 



—y 



0(e.,T)-0(e«'a) 



(B34) 



where Af is the number of lattice sites and 9{e) is the step function. The spin wave velocity Cg is given by 



1 1 



■E 



(B35) 



These expressions can be evaluated using the numerical soliton solution in 14 x 14 lattice for the one-hole problem. 
The difference of the stiffness from that at half-filling can be viewed as the renormalization due to the hole. In 
the presence of a finite density of holes, the naive generalization is to multiply the one-hole renormalization by the 
hole density. However, caution must be exercised if one attempts to apply Eqs. (B34) and Eqs. ( B35| ) to large V jt 
ratios. This is because these two expressions involve not only the ground state but also all electron-hole pair excited 
states. It is reasonable to expect that the mean field solution is a good approximation to the ground state since the 
self-consistency has been achieved. But the excited states constructed by creating electron-hole pairs in the mean-field 
spectrum may not be a good approximation of the true charge-excited states of the full Hubbard model. 

So far we have calculated the Gaussian fluctuations around the saddle point solution, i.e. the soliton. Therefore, only 
small fluctuations have been taken into account. Since there are many degenerate saddle point solutions, corresponding 
to the different locations of the solitons, wc must sum up the fluctuation contributions around all saddle point solutions. 
Furthermore, we must also include the tunneling processes between the different saddle point solutions. In terms of 
the eigenoperators of the saddle point Hamiltonian corresponding to (BIO) whose Ai*'''(i) profile has a dip located at 
the site lo, the soliton creation operator ]\ is defined through 



|Bag(*o)) - n iVacuum ) '^'^^"^ ^ |SDW) , 



(B36) 



Lower bands 



where |SDW) is the state without hole. The soliton operator has finite overlap with the bare hole operator in Eq. (|B4|), 

Qo,<To = \/^bag /I + •••• (B37) 



For a given site, uo can only take one value depending on the direction of the magnetization on the site iq. Note that 
the SU(2) rotational symmetry is restored by the fluctuating g{i). Numerically, it is straightforward to verify, in the 
14 X 14 lattice, that X^i <t ^(^cr^icr) — —1 and \ ^i^{Sf)\ — 1/2. Since all the spatial variations take place locally 
around the center of the soliton, we have to assign these quantum numbers to the soliton. Similarly, we can take 



another pair of Ai°^(i) and t^c\i) profiles with exactly the same spatial variation except that the dip is centered 



.(0) 



on a nearest neighbor site of iq. Thus, we construct another soliton in complete analogy with Eq. (B36) and define 



the soliton creation operator /' on the nearest neighbor site of ig. If we repeat the calculations determining the 



quantum numbers of the soliton, we find ^ 



-1, same as before, but J^i^i^i) changed sign. So 



the two solitons const ruc ted above have opposite quantum numbers. The interaction part of the Hamiltonian 
corresponding to Eq. (B5) has a nonzero matrix element between these two solitons: 



(Bag(jo)|<5/:|Bag(zo)) [5(^0)5+ (jo)] 12 , 



(B38) 



where t = t x [overlap integral! . This corresponds to a hopping term of the soliton in the effective action (B7), 
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Using the representation of the SU(2) matrix g{i), i.e. Eq. (p2|), we can rewrite [g{io)g~^Uo)]i2 = ^i(*o)-22 (jo) — 
Z2{'io)zi{jo)- Thus, we see that the tunnehng term between neighboring sohtons is precisely the hopping term of the 
t-J model when it is written in terms of the slave fermions and Schwinger bosons in conjugate representations for the 
two sublattices. Apart from renormalizing t to t, the only difference is that the correlation between the hole location 
and the spin magnitude has been included in the mean-field calculation of the Hubbard model. This correlation has 
been taken into account in determining the size of the soliton. 

So far, we have eliminated the high energy process es of order U. The new "elementary particle" is the soliton. The 



effective action at this level consists of Eq. (^) and ( B39 ) which are equivalent to the t-J model (^3|). The remaining 
energy scales are t (we redenote t by t) and J. In the next step, we want to construct the quasiparticle from the soliton 
and eliminate the energy scale t. Two of the present authors have recently shown how to construct a quasiparticle state 
in the Neel state and calculated the quasiparticle dispersion and spectral weight .EJ When the Neel order is destroyed 
by finite doping, we shall construct the quasiparticle state from a smoothly varying antiferromagnetic background. 
The spirit is the adiabatic approximation: The hopping of the solitons is a local and high-energy process while the 
low-energy spin excitations are described by slow variation of the spin directions. Without holes, the short wavelength 
spin fluctuations cost a lot of energy and can be eliminated. The result is the renormalization of the coefficients in 
the non-linear sigma model. In the doped case, the energy cost of short wavelength spin fluctuations near the solitons 
is compensated by the energy gain resulting from the hopping of the solitons. Therefore, short wavelength spin 
fluctuations cannot simply be eliminated as in the undoped case. 

Let \SRO) be a typical spin configuration in a short range antiferromagnetically correlated state in an arbitrary 
region of space, this configuration can be obtained from the perfectly ordered state by applying smoothly varying 
SU(2) rotations: 

\SRO)=l[g{i)\SDW). (B40) 

i 

We apply the soliton creation operator to this configuration to generate a soliton at an arbitrary site 

\to,l = 0)=fl\SRO), (B41) 

where we used the notation Z = to indicate that the number of highly frustrated spins is zero, i.e. there is no spin 
which is almost parallel to its neighboring spins. Applying on jig, I — 0) the effective Hamiltonian for the solitons, i.e. 
the t-J Hamiltonian (p3|), a new state \iQ,l = 1) is generated, 

Ht-j\io, l = 0)= -2t |^o, / = 1) + (^o, Z = 0\Ht-j\io, I = 0) |^o, I ^ 0) + ■ ■ ■ , (B42) 

where \io,l = 1) is generated by the hopping term of the t-J model. It is a superposition of four configurations in 
which the soliton has made one hop from the original site io to anyone of the four neighbors. The notation I — 1 
indicates that there is one spin, at the site io, which is frustrated with its neighbors. This is the generalization of the 
string state of one unit length discussed in Ref. ^ to the smoothly varying antiferromagnetic background. The other 
generalized string states, \io,l) for I > 1, are constructed by repeatedly applying the t-J Hamiltonian to \io,l — 1). 

Ht-j\io, 1^1)^ -2t\io, 1^0)- V3t\io, I = 2) 

+ {io, I = l|i^t„J|^o, I = 1) |^o, Z = 1) + • • • , (B43) 
Ht-j\io, I) = -V3t (Ko, + \io,l + 1» + {io, l\Ht-j\io, No, + • • • , l>2. (B44) 



Eqs. ( |B42| ) , (B43) and (B44) are the effective representation of the t~J model in the determination of the quasiparticle 



structure. Diagonalizing these equations with appropriate (io, ^l^^t- j|*Oj gives rise to a quasiparticle at the site io 
of the form, 

oo 

|io) =^u/No,0, (B45) 
;=o 

where ui are numerical coefficients which decrease exponentially as / increases. The quasiparticles located at other 
lattice sites are constructed in the same way. 

The formation of the quasiparticle further renormalizes the coefficients in the non-linear sigma model. Although it 
is not easy to calculate the renormalization accurately, it is rather easy to see that hole (or soliton) hopping inside the 
quasiparticle reduces the spin stiffness. The spin stiffness is roughly proportional to the longitudinal spin magnitude. 
The soliton hopping inside the quasiparticle reduces the average spin magnitude (S) by an amount of {2{1){S) + 1) 



20 



times the density of holes, where (l) is the average string length inside the quasiparticle, Eq. (B45). For J/t < 1, we 
have (/) > 1. 

Now, we show that the hopping amplitude of the quasiparticle constructed in the way described above contains a 
Berry phase factor. Let \i) and |j) be two quasiparticle states at two different sites of the same sublattice, i — j = x + y 
or 2x, we want to show that the quasiparticle hopping amplitude in the tight-binding approximation contains the 
appropriate phase factor, i.e., 

mt-j\j) - {t\Ht-j\i){t\j) - a{i^j){n,\n,), (B46) 

where a{r) is a real coefficient, and \nj) is the coherent spin state defined through n ■ S\n) — S\n). Since the 
quasiparticle hopping is generated by the quantum spin fluctuations, it suffices to show that the matrix element of 



the Heisenberg operator {i\Hj\j) contains the appropriate phase factor. Using Eq. (B45), we write 



(ili/jlj) (B47) 

The state |j, is a superposition of spin configurations; Each of them contains a string of length So the matrix 
element {i, l\Hj\j, I') is a sum of matrix elements of Hj between two string configurations. We need to show that the 
matrix element of Hj between any two string configurations contain the phase factor {ni\nj). In the Bethe lattice 
approximation(without winding pathL_the spin exchange operator in Hj generates quasiparticle hopping by chopping 
off the strings in by two units.t2l In all these processes, we can explicitly verify that {i,l\Hj\j,l') oc {ni\nj). 

A general example is shown in Fig. |l^. The existence of winding paths (Trugman process) leads to additional 
contributions to the matrix elements {Ll\Hj\j,l'). In general, these contributions are not exactly proportional to 
{ni\nj). An example is shown in Fig. |lj. Although not exact, the phase factor is approximately proportional to 
{fiilfij), for smoothly varying spin directions. Furthermore, the Trugman process occurs only for I + I' > 6 and 
therefore it is multiplied by a small coefficient uiuii. Thus, we have shown that the quasiparticle hopping amplitude 
acquires the Berry phase factor if the spin directions vary smoothly in space and time. The Berry phase factor can 
be expressed in terms of the SU(2) matrices g{i). Depending on which sublattice i and j belong to, the Berry phase 
factor (fiilnj) is either [g{i)g~^{j)]ii oi' [ff(*)5^(j)]22- quasiparticle hopping between different sublattices, 

i — j = 2x + y OT 3x, the appropriate phase factor of the hopping amplitude can be determined in the same way. 
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FIG. 1. Cartoon picture of the quasiparticles, represented by the smaU spheres, surfing on the fiuctuating membrane 
whose normal directions define the local directions of the magnetic moments. The north pole of the spheres represents the 
quasiparticle spin direction which must be perpendicular to the membrane. The normal vectors of the membrane form a sphere 
of unit radius. The quasiparticle traveling a closed path picks up a phase that is proportional to the area in the unit sphere 
enclosed by the normal vectors of the membrane along the quasiparticle trace. 

FIG. 2. Hartree-Fock solution of a holon in a 100-site chain with open boundary condition for U /t = 8/3. The zigzag of the 
spin density near the two ends of the chain is due to the open boundary condition and has no infiuence on the holon solution 
at the middle of the chain. 

FIG. 3. The one-electron spectrum of the Hamiltonian (^) for U /t = 8/3. For clarity the whole spectrum has shifted up by 
amount 0.148, i.e. taking fi — fi ~ U/2 = —0.148 in Hamiltonian (E6|). There are two localized states inside the Mott-Hubbard 
gap. Each localized state contains 50% spectral weight transferred from each of the two Hubbard bands. 

FIG. 4. The wave functions of the midgap localized states of the one-electron spectrum around the holon for U/t = 8/3. 
The lines are the guide to the eye. 

FIG. 5. Fitting the spin and charge densities of a holon to analytic functions for U/t — 8/3. The points are the 
numerical solution. The lines are (rii) — 1 — l/(2^c cosh^[(i — jo)/?c]) with ~ 3.209 for the charge density, and 
2(-l)'(S'f > = 0.545 tanh[(i - io)/^s] with = 2.978 for the spin density. 

FIG. 6. The magnetization profile around the kink of optimal size for U/t — 8/3. The length of the arrows is proportional 
to the spin expectation value (S^). 

FIG. 7. The coarse grained magnetization profile of a holon. The length of the arrows is proportional to the spin expectation 
value {{Sf + Sf+i)/2), and the arrows are placed at the middle of the lattice links. After coarse graining, the magnetizations 
are nonzero only around the kink and they sum to zero. Thus, the holon carries 5^ = 0. 

FIG. 8. Variational energy of the kink when an additional electron is filled into the midgap states for U/t — 8/3. The kink 
as a local deformation becomes unstable and decays into a smooth twist. The reference energy is that of the SDW state at 
half-filling of the 100-site chain with open boundary condition. The lines are the guide to the eye. 



FIG. 9. Sketch of a spinon excitation which is a smooth twist of the local magnetizations. 

FIG. 10. The profile of the spin and charge densities 2( — l)'"'"*"'" (S'f ) and 1 — (rii) of the spin bag solution for U/t — 6 in a. 
14 X 14 lattice with periodical boundary condition. Only the region near the bag is shown. The lattice spacing is the same as 
the mesh unit. 

FIG. 11. The one-electron spectrum of the spin bag solution for U/t — 6. There are two localized states inside the 
Mott-Hubbard gap. Each midgap level can be thought of as being pushed out of one of the two Hubbard bands. 

FIG. 12. The variations of the local magnetizations from the half-filled case for U/t = 6. The length of arrows is proportional 
to S{S^). The length of the center arrow has be reduced by a factor of 10 for clarity. The changes are localized around the 
center of the spin bag, and they sum to ±1/2. The sign depends on which sublattice the center of the spin bag resides. Thus, 
the spin bag carries — ±1/2, in contrast to the holon in ID which carries — 0. Note that the variation of the staggered 
magnetization has the same meaning as the coarse grained ones shown in Fig. ^ since the coarse graining at half-filling over 
each plaquette yields zero everywhere. 
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FIG. 13. A general example of the processes generating quasiparticle hopping with no winding path. Starting from 
configuration A, we remove a hole either at the site i or j, marked by the circles, to generate B or B'. The hole in B hops three 

times to generate C which contains a string of length 3 marked by the dashed line. C' is generated from B' by one hop of the 
hole. C and C' belong to \i, I = 3) and \ j, I' = 1) respectively. Hj contains the spin exchange operator that can exchange the 
first two spins along the string in C. The resulting configuration after exchange is shown in D together with C' . The matrix 
element {i,l = S\Hj\j,l' = 1) is proportional to {ni\nj). 

FIG. 14. An example of the Trugman process generating quasiparticle hopping. The configurations B and B' are generated 
by removing a hole in A either at the site a or e, marked by the circles. C is a configuration in |o, 1 = 6). Hj can exchange 

spins e and b in configuration C, generating D. The contribution of this example to {a, I = 6\Hj\e,l' = 0) is given by the 
overlap between D an B' , and is proportional to {na\nc) {nc\ne) — (na\ne) in the smooth antiferromagnetic background. 
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